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Abstract 

We consider a ferrofluid cylinder, that is rotating with constant rotation frequency ft = £le z as 
a rigid body. A homogeneous magnetic field Ho = Hqg x is applied perpendicular to the cylinder 
axis e z . This causes a nonequilibrium situation. Therein the magnetization M and the internal 
magnetic field H are constant in time and homogeneous within the ferrofluid. According to the 
Maxwell equations they are related to each other via H = Ho — M/2. However, H and M are 
not parallel to each other and their directions differ from that of the applied field Ho- We have 
analyzed several different theoretical models that provide equations for the magnetization in such 
a situation. The magnetization M is determined for each model as a function of £1 and Hq in a 
wide range of frequencies and fields. Comparisons are made of the different model results and the 
differences in particular of the predictions for the perpendicular components H y = —M y /2 of the 
fields are analyzed. 
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I. INTRODUCTION 

There are several theoretical equations for the dynamics of the magnetization M(r, t) 



experimentally accessible. We consider a ferrofTuid cylinder of radius R of sufficiently large 
length to be approximated as infinite in a homogeneous applied field H = H Q e x in x- 
direction. The ferrofluid cylinder is enforced via its walls to rotate as a rigid-body around 
its long axis with constant rotation frequency O = Qe z being oriented perpendicular to Ho. 
The flowfield is thus u(r) = f2 x r = Vtre.^ where is the unit vector in azimuthal direction. 
In such a situation all aforementioned models allow for a spatially and temporally constant 
nonequilibrium magnetization M that is rotated out of the directions of H and H by the 
flow. The Maxwell equations demand that the fields H and M within the ferrofluid are 
related to each other via 



as indicated schematically in Fig. ^ and that the magnetic field outside the ferrofluid cylinder 




H = H - - 

° 2 



M 



(1.1) 




(1.2) 



is a superposition of the applied field Hq and the dipolar contribution from M. 




FIG. 1: Schematic plot of relevant vectors. 
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II. MAGNETIZATION EQUATIONS 

The model equations that we compare here imply a relaxational dynamics either of M 
towards the equilibrium magnetization 

M eg (H) = M^1 H = x{H)H (2.1) 

or of the "local equilibrium" or "effective" field 

H e// (M) = eq ^ J M = F(M)M (2.2) 

towards the internal field H. The equilibrium magnetization M eq (H) referring to the func- 
tional relation between internal field H and magnetization in the case of Q = is a thermo- 
dynamic material property of the ferrofluid. The effective field H e ff lies parallel to M and 
can be seen as the inverse of the defining requirement 

M = M e? (H e// ) . (2.3) 

In equilibrium, Q = 0, one has H e ff — H and M = M eg . 
We consider here the relations 

- (M - M eq ) (2.4) 

r 

i (M - M eq ) + x (M x H) (2.5) 
lH (H e// - H) + f^M x (M x H) (2.6) 

S'01 [4(:_ " x H e// = i (H e// - H) + ^H e// x (M x H) (2.7) 



Debye: 


ft x M 


S'72 Q: 


x M 


fk y : 


n x m 



ML 0: fi x M = f(H e// - H) (2.8) 

resulting for the rotating cylinder from the above 5 models. In ML we use the weak field 
variant of ref. Q]. These equations have to be solved numerically in combination with the 
Maxwell equation Ql.lj) . 

As an aside we mention that the above equations can be written in the common form 

M x (n + a 3 M x H ) = «i(H - a 2 M) (2.9) 



with coefficients: 
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III. RESULTS 

In order to make the comparison of the theoretical results easier we replace the equilibrium 
magnetization M eq (H) by the Langevin expression M eq (H) = M sat C (3xoH/M sat ) with the 
initial susceptibility xo — x(H = 0). We use Xo = 1-09 and M sat = 18149A/m for the 
saturation magnetization which is appropriate for the ferroffuid APG 933 of FERROTEC. 
The resulting curve is shown in Fig. El Furthermore, we replace the relaxation time t(H) 
by Tg = 6 ■ lCT 4 s. For ( ~ |$r] we use the values rj = 0.5Pa ■ s and $ = 0.041 and for ^ H 
we use = Xo/tb fl- For the parameter ^ of ML \^ we investigate two different choices: 
Either the low-field variant, £ = Xo/ t b, & s in FK that is denoted here by ML(F). Or the 
variant £ = 1/[F(M)tb] as in S'Ol that is denoted here by ML(S). 

Especially the perpendicular component H y = —\M y of the magnetic field is suited for a 
comparison of the different models with each other and with experiments. Before doing the 
former we should like to draw the attention to the frequency behavior of M y (H ,Q). We 
mentioned already that M y vanishes for zero vorticity, Q = 0. Furthermore, one finds that 
M y as well as M x vanishes also in the limit Q — ► oo. And since one can rewrite the solution 
of eq. (|2.9j) in the form M y = ai _^ M 2 jj^ one sees that M y (Q,) has a maximum as a function 
of Q as in Fig. |3J There we show H y versus Q. 

The differences in the results for the different models are easily captured by comparing 
their predictions for the maximum values of \H y \, the locations of these maxima at Q max , 
and the initial slopes at — > 0, each as a function of applied field H . This is done in 
Fig.H 

The maximal values of \H y \ of Debye and S'72 are the same while their locations, Q max , 
differ. The models S'Ol, FK, and ML formulated in terms of the effective field also share a 
common maximal value of \H y \ being larger than that of Debye and S'72 while the location, 
n max , differ partly substantially. Hence the magnetic torque, M x H, entering into S'72, 
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FK, and S'Ol only shifts the frequency Q max . It remains to be seen whether experiments 
can be performed with sufficient accuracy to discriminate between the different theoretical 
predictions. 
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FIG. 2: Equilibrium magnetization M eq (H) used as input into the models compared here. 




FIG. 3: Comparison of the predictions of the different theoretical models for the transverse internal 
field H y versus rotation frequency VL. 
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FIG. 4: (a)Frequency Q max leading to maximal transverse field, (b) largest transverse field, and 
(c) initial slope at SI — > 0. 



